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We study the dynamics of a Bose-Einstein condensate in 
a double-well potential in the two-mode approximation. The 
dissipation of energy from the condensate is described by the 
coupling to a thermal reservoir of non-condensate modes. As a 
consequence of the coupling the self-locked population imbal- 
ance in the macroscopic quantum self-trapping decays away. 
We show that a coherent state predicted by spontaneous sym- 
metry breaking is not robust and decoheres rapidly into a sta- 
tistical mixture due to the interactions between condensate 
and non-condensate atoms. However, via stochastic simula- 
tions we find that with a sufficiently fast measurement rate 
of the relative phase between the two wells the matter wave 
coherence is established even in the presence of the decoher- 
ence. 
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Bose-Einstein condensates (BEC's) exhibit a macro- 
scopic quantum coherence which in thermal atomic en- 
sembles is absent [Q. In the conventional reasoning the 
BEC is assigned a macroscopic wave function with an 
arbitrary but fixed phase. The selection of a phase 
implies the spontaneous breaking of the gauge symme- 
try. The atom- atom interactions in finite-sized BEC's 
affect coherence properties. The relative phase of BEC's 
undergoes diffusion or collapses due to the condensate 
self- interactions jHj|. The interactions between conden- 
sate and non-condensate atoms create decoherence [|J. 
Modelling decoherence by fully including the quantum 
effects requires sophisticated theoretical studies which 
non-trivially include non-condensate atoms. In the ex- 
periments on BEC's of dilute alkali atomic gases ||] 
trapped atoms are evaporatively cooled and continuously 
exchange particles with their environment. Thus, stan- 
dard approaches of quantum optics for open systems in- 
volving master equations and heat reservoirs seem espe- 
cially natural for treating atomic BEC's [p|-|io[. 

In this paper we study the evolution of the master 
equation for a BEC in a double-well potential in a two- 
mode approximation using previously derived models . 
Macroscopic quantum coherence of BEC's results in co- 
herent quantum tunneling of atoms between the two 
modes or the "two BEC's", which is analogous to the 
coherent tunneling of Cooper pairs in a Josephson junc- 
tion 1 1 - 1 5 1 . According to the Josephson effect , the atom 
numbers of the BEC's oscillate even if the number of 
atoms in each condensate is initially equal. Even BEC's 
with a well-defined number of atoms, and with no phase 
information, could exhibit oscillations in particular mea- 



surement processes on atoms |L6| or on photons [ pL 7|Jl8| . 

One interesting feature of the coherent quantum tun- 
neling between two BEC's is that due to the nonlinearity 
arising from atom-atom interactions, oscillations are ex- 
pected to be suppressed when the population difference 
exceeds a critical value in a process known as macro- 

II. We show 



scopic quantum self-trapping (MQST) |l£ 
that in the presence of collisions between the condensate 
and non-condensate atoms MQST decays away, i.e., the 
atom numbers of the two BEC's become balanced. 

Anglin derived the master equation for a trapped 
BEC by considering a special model of a BEC confined in 
a deep but narrow spherical square-well potential, tuned 
so as to posses exactly one single-particle bound state. 
In this model the thermal reservoir of non-condensate 
atoms consists of a continuum of unbound modes ob- 
tained by the scattering solutions of the potential well. 
The binding energy Eb of the BEC mode is assumed to 
be large compared to the thermal energy of the noncon- 
densed atoms. Then the fugacity satisfies z — <C 1, 
where p is the chemical potential and (3 = l/(ksT), even 
in the presence of a BEC. The condensation occurs be- 
cause of the depth of the attractive potential. The small 
fugacity allows the derivation of the master equation in 
the Markov and Born approximations. By expanding in 
terms of the small parameters z and a/d, where a is the 
s-wave scattering length and d is the length scale of the 
BEC, the reduced density operator for the BEC satisfies 
the following equation of motion: 

p = i/h[p, H S ] + CiD[tfb]p + C 2 V[tf]p 

+C 2 exp {(3(hA - p)/2}V[b]p + 0(z 2 a 3 /d 3 ) . 



(1) 



Here we have defined 



T>[c]p = cpc< — l/2(c^cp + pe'e) , 



(2) 



and A ~ 2kN — E\, is expressed in terms of the strength 
of the self-interaction energy k. The system Hamiltonian 
for the BEC is denoted by Hs- For simplicity, the inter- 
actions between different non-condensate atoms are esti- 
mated by the Boltzmann scattering rate 7 = anhk/m, 
where a is the scattering cross-section and n is the den- 
sity of the gas. The parameters C\ and C2 may then be 
estimated by C\ ~ 7 and C2/C1 ~ z/((3Eb). 

In a single mode approximation the processes in which 
a BEC atom collides with a non-condensate atom and 
produces two BEC particles, or vice versa, do not con- 
serve the energy and they are absent in Eq. (0). The 
term proportional to C\ describes elastic two-body colli- 
sions between condensate and non-condensate atoms and 
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it induces phase decoherence |T|] . Inelastic collisions [the 
terms proportional to C2 in Eq. (|l|)] introduce amplitude 
decoherence [Gl|. With the present approximations the 
amplitude decoherence is dramatically reduced compared 
to the phase damping. The central assumption is a small 
fugacity indicating /3Eb ^> 1. In the scattering processes 
the energy must be conserved to leading order, so that 
depletion and growth of the BEC involve non-condensate 
atoms with high enough energy to balance the large bind- 
ing energy of the trapped state. If the BEC is described 
by a multiple number of modes the amplitude de- 
coherence does not need to be small compared to the 
phase decoherence. The necessary condition for the va- 
lidity of the single-mode approximation in a harmonic 
trap is that the self-energy of the atom-atom interac- 
tions of the BEC should not dominate over the mode 
energy spacing. This assumption clearly breaks down 
in the Thomas-Fermi limit in which the kinetic energy is 
negligible compared to the self-interaction energy indicat- 
ing (157Va/0 1/5 > 1, where I = (Ti/muj) 1 / 2 is the length 
scale of the harmonic oscillator. Jaksch et al. || have cal- 
culated the intensity and the amplitude fluctuations of a 
BEC. They have evaluated the coefficients of the master 
equation in the Thomas-Fermi limit and have obtained 
stronger amplitude decoherence than phase decoherence. 
Nevertheless, the calculations of the BEC fluctuations 
are still performed in the one-mode approximation. 

We consider Anglin's model for the master equation 
in the studies of the coherent quantum tunneling of 
a BEC in a double-well potential. To obtain the system 
Hamiltonian for the BEC we approximate the total field 
operator by the two lowest quantum modes and write 
ip(r) ~ </>fc(r)6 + c (r)c, where 4>b and </> c are the local 
mode solutions of the individual wells with small spatial 
overlap and the corresponding annihilation operators are 
b and c. This requires that the BEC self-interaction en- 
ergy should not dominate over the energy separation of 
these modes from the higher excited modes. The Hamil- 
tonian in the two-mode approximation reads |l3f 

H s /h = £tfb + Q(& f c + c f 6) + n[(tfj 2 b 2 + [c^fc 2 ] . (3) 

Here £ is the energy difference between the modes. The 
tunneling between the two wells is described by Q, which 
is proportional to the overlap of the spatial mode function 
of the opposite wells. The short-ranged two-body inter- 
action strength is obtained from k = 2irah/m J |</>f,(r)| 4 , 
where we have assumed J \4>b{r)\ 4 = f |0 c (r)| 4 . 

It is useful to describe the dynamics of Eq. (||) in terms 
of the atomic coherent states p0[ , which exhibit coher- 
ence, but conserve the total number of atoms: 
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where |j, m) is an eigenstate of the angular momentum 
operators J z and J 2 with eigenvalues m and Here 
r = e~ 1 ^ tan (0/2). We define the angular momentum 



operators in terms of the BEC operators in the usual 
way: J x = (tfc + c)b)/2, J y = (tfc- c*b)/{2i), and J z = 
(b^b — c^c)/2. Then j — N/2 is a constant of motion and 
the atomic coherent states may be expressed in terms of 
the number states of the BEC's 
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where / 



m- 



-N/2, and B and C are the 'coherent ampli- 



tudes' of the two BEC's: (tfb) = \B\ 2 , (c+c) = |C| 2 , and 
(tfc) = \B\\C\e lip . The relative phase between the two 
wells is ip. Equation ([5J) clearly shows how the atomic co- 
herent states are projections of the coherent states \B,C) 
onto the basis sets of a fixed total number of atoms N. 

We study the evolution of atomic coherent states in 
the presence of decoherence in both wells b and c. As 
already noted, with the present approximations the am- 
plitude decoherence is much weaker than the phase de- 
coherence. It is also easy to verify numerically that the 
coherent states are much less robust for the phase deco- 
herence than for the amplitude decoherence, even if the 
magnitude of the phase and the amplitude damping is 
equal. This is also well known in quantum optics |l9| ]. 
Hence, the effects of the amplitude decoherence are neg- 
ligible compared to the phase decoherence. 

The dynamics of the master equation is studied in 
terms of stochastic trajectories of state vectors pi| . The 
master equation is unraveled by Monte-Carlo evolution 
of wave functions. 

In Fig. pi) we have plotted the expectation value of the 
number of atoms in the well 6, Nb(t). The initial state is 
the atomic coherent state with the relative phase ip = ir/2 
between the two wells and the expectation values for the 
atom numbers Nf, = N c = 50. We have set Nn/Sl = 0.5 
and £ = 0. The solid line is the result without deco- 
herence 7 = 0. The oscillations are damped due to the 
collapse of the macroscopic coherence [Oi. The dashed 
line has 7 = 0.1k and the dashed-dotted line 7 = 0.4k. 
The decoherence clearly increases the damping of the os- 
cillations. Although the model used is very simplified, we 
can make rough estimates for the parameters. For the ef- 
fective mode volume 1/ J |0(r)| 4 = 10~ 9 cm 3 , a = 5 nm, 
for 23 Na, and for the temperature T = 100 nK, 7 = 0.1k 
corresponds to the density of the non-condensate atoms 
n ~ 10 9 cm~ 3 and the fugacity z ~ 10~ 3 . In Fig. |l|b) we 
have plotted Tr(p 2 ) for the same run of simulations. We 
see that a pure state predicted by spontaneous symmetry 
breaking is not robust and decoheres rapidly into a statis- 
tical mixture due to the interactions between condensate 
and non-condensate atoms. 

If the nonlinearity is large compared to the tunnel- 
ing frequency and the population imbalance exceeds a 
critical value, the oscillations of the atom numbers are 
suppressed |l3| , |l4| ]. A large number of atoms remains 
"locked" in one of the wells. In Fig. 0a) we have plot- 
ted Nb(t) obtained by integrating Eq. (0) (the solid line) 
and the solution of the master equation in the presence 
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of the decoherence 7 = 0.2k (the dashed line). In this 
case Nk/VL = 4.5, £ = 0.00517, and the initial state is 
the atomic coherent state with the expectation values 
Nb = 4N C = 80 and ip — 0. Due to the interactions 
between condensate and non-condensate atoms MQST 
vanishes and the atom population becomes balanced. 

Next, we include the effect of measurements into the 
calculations. We assume that the number of atoms is 
nondestructively measured in one of the two wells. The 
effect of the measurement is included in quantum trajec- 
tory simulations by averaging over the dissipation chan- 
nels corresponding to the interactions between conden- 
sate and non-condensate atoms, but at the same time by 
considering the measurement of the number of atoms in 
one of the wells to be a single realization of a stochastic 
trajectory. We consider a particular situation in which 
the Josephson dynamics is nondestructively measured by 
shining a coherent light beam through one of the BEC's. 

We assume that the incoming light field with a large 
detuning from the atomic resonance is scattered from the 
well b. For instance, if the shape of the gas is flat and the 
light is shined through a thin dimension, the dipole shifts 
are small and the sample can be considered optically thin 
p2fl . A BEC atom scatters back to the BEC via coher- 
ent spontaneous scattering, stimulated by a large num- 
ber of atoms in the BEC. Coherently scattered photons 
are emitted into a narrow cone in the forward direction. 
The amplitude of the scattered field has the dependence 
|Eg| oc £d eg /(hA) b% on the detuning A, the dipole ma- 
trix element d eg , and the amplitude of the incoming field 
£ p2| |. The direct counting of spontaneously emitted 
photons can be simulated in terms of quantum trajecto- 
ries 21 1, in which the stochastic quantum "jumps" cor- 
respond to the detections of photons. The procedure is 
similar to Ref . [^3| . The detection rate of the scattered 
photons in the present case is T((b^b) 2 ) oc |Eg • E^|. 

If the number of atoms in a BEC is not large, the 
scattering between the condensate and non-condensate 
modes is not negligible. This introduces amplitude de- 
coherence similar to the amplitude decoherence due to 
the atomic collisions in Eq. (Q). If we require that the 
two-mode approximation describes accurately the tun- 
neling dynamics, for small harmonic traps the amplitude 
decoherence due to the light scattering may not be neg- 
ligible. For large traps the two-mode approximation can 
be accurate even for large atom numbers as the BEC 
self-interaction energy k oc l~ 3 and the trap frequency 
lu oc l~ 2 . Nevertheless, as a first approximation we ig- 
nore the decoherence due to the light scattering. 

The density matrix of a BEC may be reconstructed by 
the nondestructive measurements of the number of atoms 
in one of the wells. The procedure is similar to Ref. 
except that the density matrix has now time dynamics 
determined by the Hamiltonian (J^), the dissipation, and 
the back-action of the measurements. Following the nota- 
tion of Ref. {H, at the time t we have p(t) = U{t)^pU{t), 
where U (t) is in this case the time propagator. Then the 



probability of the measurement result of m atoms in the 
well b at the time t is given by P m (t) = (m\p{t)\m) . By 
inverting this expression the density matrix can be re- 
constructed |p4) . 

The off-diagonal long range order (ODLRO) between 
the two wells may be described by the visibility of the in- 
terference [3 (It]] . To emphasize the effect of decoherence 
on p we ignore the oscillating dynamics of H$ (includ- 
ing the collapses and revivals) by propagating back the 
system dynamics. In accordance with Ref. |T^] we define 



- = ?-Tr[e lHst/h pe~ tHst/n b f (0)c{0)] . 



(6) 



For a coherent state we have j3 = 1, and ip is the relative 
phase between the two wells. However, for a number state 
there is no phase information and (3 — 0. If the BEC's 
have unequal atom numbers, the maximum visibility is 
reduced from one to /3 m ax = 2^/NbNc/N. Hence, it is 
useful to define the relative visibility by f3 r = /3//3 max . 

We simulate the dynamics of the dissipation and the 
measurements by repeating single realizations of quan- 
tum trajectories. In the first realization we save the 
stochastic times of the photon detections. In every sub- 
sequent run of the trajectory the times of the photon 
detections are forced to be the same as in the first run. 
Although the photon detection times after the first tra- 
jectory are deterministic, the collision times between con- 
densate and non-condensate atoms corresponding to the 
dissipation channels are stochastic in every run. Aver- 
aging over all the trajectories allows us to consider the 
photon measurements to be a "single realization" of the 
quantum trajectory even though the atomic collisions are 
at the same time ensemble averages corresponding to the 
density matrix evolution. 

We consider a situation that the two BEC's are ini- 
tially in pure number states with iV& = 52 and N c = 48. 
We have set JVk/U = 0.25, £ = 0.005rj, and the photon 
scattering rate T = 0.8k. In Fig. ||a) we have plot- 
ted (3 r (t). The solid line is the result without decoher- 
ence 7 = 0. The dashed line has 7 = 0.05k and the 
dashed-dotted line 7 = 1.8k. In the beginning for the 
number state (3 r = 0, but (3 r — > I rapidly even in the 
presence of weak decoherence. As a consequence of the 
decoherence ODLRO starts decreasing, but the measure- 
ments of spontaneously scattered photons establish the 
macroscopic coherence, even though the BEC's are ini- 
tially in pure number states. In Fig. ||b) we have plotted 
Tr(p 2 ) for the same run. We see that Tr(p 2 ) remains close 
to one and the state is reasonably pure due to the fast 
measurement rate even in the presence of decoherence if 
7 = 0.0625r. In the case of stronger decoherence with 
7 = 2.25r the state evolves into a statistical mixture. 

In conclusion, we have shown that as a consequence of 
the interactions between condensate and non-condensate 
atoms MQST decays away. Due to the interactions a 
BEC does not remain in a pure state with a well-defined 
relative phase. However, the coherence properties can 
be established via the measurement process even in the 
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presence of decoherence. In particular, nondestructive [14 

detections allow the measurements of phase dynamics. [15 
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FIG. 1. The expectation value of the number of atoms in [18 
the well b in a) as a function of time. Initial state is the atomic [19 
coherent state with Nt — N c = 50 and the relative phase 
ip — tt/2. The solid line is the result without decoherence [20 
7 = 0. The dashed line has 7 = 0.1k and the dashed-dotted [21 
line 7 = 0.4k. In b) Tr(p 2 ) for the same run. 

[22 

FIG. 2. The expectation value of the number of atoms in [23 
the well b in the case of large nonlinearity. Initial state is the [24 
atomic coherent state with Nt ~ 4iV c = 80 and ip = 0. The 
solid line is the result without decoherence and describes the 
macroscopic quantum self-trapping. The oscillations undergo 
collapses and revivals. The dashed line shows how the atom 
population becomes balanced in the presence of decoherence. 
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